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Abstract We demonstrate that the multiple gluon emission phase space in the dipole 
cascade model has a strong linear correlation with the number of gluons emitted. The 
number of gluons per unit available phase space at a certain resolution scale is found to 
be remarkably independent of the cms energy and global event properties like thrust, and 
even changes in the ordering variable or resolution scale. We show that the distribution of 
sizes of gluon-gluon dipoles in a parton cascade has stability properties which are sufficient 
to account for those of the phase space variable. We observe that certain more abstract 
entities, defined in the context of hadronisation and related to the gluon emission phase 
space, share those properties of colour dipoles and name them Generalised Dipoles. We 
also present an analytical model to qualitatively describe our findings. 
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1 Introduction 



The phase space for multiple bremsstrahlung emission of photons in QED is entirely 
given by the properties of the original current. The reason is that the photon-quanta 
are uncharged and therefore the currents are not changed because of emission of photons 
(besides the recoil effects from hard emissions). 

On the other hand the QCD field quanta, the gluons, are charged. Already the emission 
of a first QCD gluon in e + e~-annihilation means that the original qq dipole is changed. 
It so happens that the change is (to a very good approximation) from one dipole to two 
independent dipoles, one between the quark (q) and the first gluon and the other between 
the gluon and the anti-quark (q). The initial dipole is at rest in the total cms whereas the 
two "new" dipoles move away from each other. The combined phase space for emitting 
from either one or the other of these dipoles is found to be larger. In rapidity space the 
increase can be described as an extra region of a size corresponding to the logarithm of 
the squared transverse momentum of the first emitted gluon. 

In the Dipole Cascade Model |IJ, as implemented in the event generator Ariadne 0, 
a second emission leads to three independent dipoles and so on. The ordering variable 
is an invariantly defined transverse momentum of the emitted gluons. The dipole masses 
decrease quickly with successive emissions. 

Herwig ||, and Pythia subdivide the emerging dipole angular regions into 
independent cones for each emitting parton. One way to state the coherence conditions is 
that they do not permit double-counting in the emission process. Herwig and Pythia 
implement the coherence conditions by means of angular ordering to avoid such overlaps. 
The ordering variable for HERWIG is just the angle (or the rapidity variable) occurring 
in the coherence conditions. Pythia uses the "virtuality" along the emission lines as an 
ordering variable and only afterwards introduces an angular ordering condition. 

In this paper we will mainly describe the emerging features in terms of the notions of 
the Dipole Cascade Model (Ariadne) and only perform a cursory comparison to similar 
results from Pythia. Our aim is to investigate a set of distributions stemming from the 
perturbative parton cascades. We will show that the cascades result in a local structure, 
corresponding to sets of independent entities, which we are going to call Generalised 
Dipoles (GDs). The GDs are linked together along the colour lines of the QCD force 
field. They have a common distribution in a generalised rapidity range A, and in the 
(local) transverse momentum. This structure is independent of the total cms energy, the 
global event variables like thrust and the number of hard emissions. Further the GDs 
show a surprisingly small scale dependence, i.e. almost the same distributions occur inside 
a wide range of the ordering variable. 

We will concentrate upon e + e~-annihilation events, where the gluon cascade is pro- 
duced through the bremsstrahlung from an originally produced qq pair (although we 
expect that the same structures will emerge also in other dynamical situations). Such a 
multi-gluon state is conveniently described by means of a four-vector valued function, the 
directrix ||. The directrix is obtained by laying out the energy-momentum vectors, 
kj, of the emitted partons in colour-order, starting with the q and ending with the q 
(we will neglect the production of "new" qq-pairs through the gluon splitting process). 
As we will use massless partons, the directrix is a curve with a tangent that is lightlike 
everywhere. 

As mentioned earlier the phase space available for gluon emission at a certain value of 
the ordering variable k± depends on the partons present and their colour order. It could, 
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therefore, be regarded as a functional of the directrix A^. In section 0, we will discuss 
a measure of this phase space, to be called £, for multiple gluon emission in perturba- 
tive QCD. An infrared stable generalisation of this phase space, involving a "resolution 
parameter" m leads to the concept of the A- measure @. 

The A- measure was introduced many years ago 0, as a generalisation of the rapidity 
range, to describe the phase space for hadronisation. At the same time |J , a curve similar 
in appearance to a set of connected hyperbolae along the directrix was introduced, the 
^-curve. Its length is related to A in the same way as the "ordinary" rapidity variable 
measures the length of a hyperbola spanned between two lightlike vectors. 

The tangent vector, to be called qx, at a point on the ^-curve that reaches out to a 
point on the directrix quickly reaches a constant length, q\ = rriQ, just like the constant 
distance in Minkowski metric between the lightcones and an ordinary hyperbola. The 
parameter mo is a resolution parameter for the properties of the directrix. In the context of 
hadronisation, it is fixed by the properties of the spectrum of hadrons produced. Thought 
of as a generalisation of £, the A-measure has a resolution parameter which can be directly 
related to the ordering variable in the perturbative cascade. When the tangent vector qx 
follows the directrix it will span a surface between the directrix and the A'-curve. The 
A-measure is also proportional to the area of this surface. 

In the dipole model, the available phase space for gluon emission, £, in the Leading Log 
Approximation (LLA) and in the Modified Leading Log Approximation (MLLA) 
schemes is a sum of a series of terms , (A£)j, one for each pair of consecutive vectors along 
the directrix ||. These terms are contributions to the phase space from individual dipoles 
in the Dipole Cascade Model. We will show that the A-measure can be similarly subdivided 
into a number of parts (AA)j, one for each vector along the directrix, and that they are 
strongly related to the contribution from the dipoles, {A£)j. The process of subdividing A 
corresponds to a partitioning of the surface mentioned above into " plaquettes" that are 
enclosed by one "initial" qx and one "final" qx vector together with a hyperbolic segment 
of the A'-curve and the parton energy-momentum vector in between, cf. Fig [I] in section 
|2j. It is the regularity in the (AA) distribution that will lead to the definition of the GDs. 
For the surface mentioned above this regularity implies a simple general structure, which 
shows only a slowly varying dependence on the ordering variable. 

The directrix also plays a major role in the description of the state of the massless rel- 
ativistic string, which is used as a model for the QCD force fields in String Fragmentation 
||. The surface spanned by the string during one period is a minimal surface. This has 
two consequences. On the one hand the surface is completely determined by its boundary 
curve. On the other hand the surface is stable against small deformations (the model is 
infrared stable). In the motion of a massless relativistic string, a wave moves across the 
string surface bouncing at the endpoints defined by the orbits of the q and the q. This 
means that the internal excitations (which in the Lund Model are identified as the gluons) 
will reach and affect the endpoints in turn, i.e. in the colour order of their emission. The 
corresponding orbit of the q endpoint is the directrix. 

It is interesting to note that the X-curve (or rather a close relative, called the V- 
curve in [|8]]) and the A-measure play a major role in the String Fragmentation process 
for multigluon string states. In the process that we devised in |8[], the final state hadron 
energy-momenta, laid out in rank order, constitute a curve, the X-curve. The X-curve 

3 In this paper, we will use the word plaquette to denote this area even though one part of the boundary 
curve is a part of a hyperbola. The plaquettes defined in this paper are not related to the plaquettes 
used in Lattice QCD. 
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follows the directrix at a typical stochastically fluctuating distance and the area in between 
the two curves could be thought of as the area in the Lund Model Area law ||. It turns 
out that the average X-curve is just the V curve with a length given by A. The generalised 
dipoles to be defined in this note also play an interesting role in the fragmentation process. 
Just as they emerge as independent entities from the partonic cascades, they also fragment 
essentially independently of each other into the final state hadrons, in a fragmentation 
scheme along the directrix. 

The A-measure has been used in many different contexts in investigations over the 
years. It was quickly found, that the <Y-curve shows properties of a fractal character || 
and the multi-fractal dimensions were shown to be identical to the so-called anomalous 
dimensions of QCD |L0] . 



In section |2|, we recapitulate some features of the gluon emission phase space ^-measure 
and the A-measure with the aim of pointing out similarities and differences between the two 
quantities^. We also introduce the quantities AA and present a geometrical interpretation 
for them. In section || we describe our findings and define the Generalised Dipoles. In 
section [|, we present a model to understand the distributions we have obtained. At the 
same time we will consider some earlier results and the theoretical analysis. Finally, in 
section |5| we make some concluding remarks. 

2 The A-Measure and The Perturbative Cascades 

2.1 The phase space in a dipole cascade and the ^-Measure 

The probability for the emission of gluons in the Dipole Model can be described in terms 
of an inclusive density of gluons (cf. [T|]): 

dk 2 

dn = a eS dy—± x (£ p ) (1) 

with a eS the effective running coupling and (J2 P ) t ne polarisation sum, i.e. the coupling 
between the spin of the emitters and the gluon. The coherence conditions in this case 
are identical to energy-momentum conservation, i.e. the cms energy of the gluon cannot 
exceed half the total cms energy 

fc±cosh(y)<^ (2) 

A convenient approximation, corresponding to the Leading Log Approximation, is \y\ < 
|ln(p-) so that the total rapidity range is Ay ~ ln(p-). As we mentioned above the 
emission of a gluon will change the current. But it is an immense simplification that this 
change is simply a change from one dipole to two independent dipoles (to a very good 
approximation) [[?], [TTH . Thus the density for the emission of two gluons is factorisable, so 
that we obtain 

dn(qgig 2 q) ~ dn(qgiq) [dn(qg 2 gi) + dn(gig 2 q)] (3) 



4 In the past, like in the references listed in this paper, the gluon emission phase space has often been 
called the A-measure. But since both this phase space and the phase space for hadronisation are relevant 
for this work, we make the distinction here and give them different names. 
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(with a negative correction that formally is of the order of 1/N% for a finite number 
of colours, but reduced further in practice due to kinematics). Eq. (||) is valid if the 
transverse momenta are ordered so that k±i > k±_ 2 or else the two gluons are exchanged. 
While the original (qq) dipole is at rest in the cms the two "new" dipoles move apart. 
Each of them will decay in its own restframe with a phase space given by Eq. (0) and this 
means that the combined rapidity region for emitting the second gluon with transverse 
momentum k i is 



(A£) 01 + (A£) 12 = ln(^ +ln ^) (4) 



The expression S01 ^ 12 is the conventionally defined invariant squared transverse momentum 
of the first gluon, k\^ Similarly, for an emission at k\ from a system of n gluons, we can 
write the phase space as: 

t = E(^W = Eln(^f i )=lnn( £ p ti ) (5) 

i=0 i=0 K ± i=0 ft -L 

At this point, it is necessary to differ between two definitions of the transverse mo- 
mentum. The one we have used up to now is the ordering variable k±, i.e. the value used 
for the emission of a gluon and defined by the partitioning of a dipole of mass M into 
a pair of dipoles with masses (mi, 777.2) so that k\ = mfm 2 ,/ 'M 2 . There is another local 
variable, to be called kt, which is defined along the directrix so that we have 



2 ^ {2kj_ikj){2kjkj + i) _ Sj—ijSjj+i 

Hl (2kj-ikj) + (2kjk j+1 ) + (2k j - 1 k j+1 ) ~~ s^-uj+i 



in terms of the energy-momentum vectors of the (colour) adjacent gluons. 

We note that the two values of k± and kt coincide for the last emitted gluon. But in 
general the transverse momentum k t of a gluon may be much smaller than the k± value 
at which it was emitted, because of the recoils from subsequent emissions. In fact, in 
this way k t can even become smaller than the resolution scale at which the state is being 
sampled. But the "emission" k± and the concepts of the "first emitted" gluon and the 
"second emitted" gluon etc are necessities of the strategies used to simulate this quantum 
mechanical process, and not properties of the final state of partons. When one asks 
what the probability is for the production of a state with a certain number of partons, 
where colour connected partons with transverse momenta smaller than a certain value 
would not be considered as resolved emissions, the only relevant transverse momentum 
is the kt, which is calculated using vectors in the final state. Therefore, it is undesirable 
that this variable is pushed to values smaller than the resolution parameter. We have 
used two different methods to avoid such situations. The first is to veto such emissions 
along the cascade[] and the second to change the ensuing directrix vectors so that three 
neighbouring lightlike vectors are combined into two. All the results of this investigation 
are independent of the method that is used. Furthermore, such a procedure is essential 



for the inclusion of the second order matrix element corrections into Ariadne 12 . 



5 This feature will be made available in the next version of Ariadne 



4 



2.2 The A-measure, its properties and its connection with the 
^-measure 

When a string without gluonic excitations fragments under the usual Lund Model as- 
sumptions, hadrons are formed by breaking of the string field at "vertices" which, on the 
average, lie on a typical hyperbola parametrised by a squared proper time from the origin, 
usually called T or ttiq. The rapidity range available along this hyperbola is given by 
In -V 

In the Lund Model interpretation for a string with a single gluon excitation, there will 
be two parts of the string; one string piece spanned between the q and the g and one 
between the g and the q. Therefore there will be two hyperbolic angular ranges (Ay)oi 
and (Ay) 12 for hadron formation. The sum of these ranges is 



(Ay)oi + (Ay) 12 = In(i^) + In(^) = In(^) + In ) (7) 



We draw attention to the factor of 4 in Eq. fl7|), and its absence in Eq. (§). In the 
Lund Model, when the partons move apart stretching a string-like field between them, 
they lose energy-momentum to the string field. A gluon is attached to two string pieces 
while a quark or an anti-quark is attached to one. Therefore a gluon will lose energy- 
momentum twice as quickly as a quark or an anti-quark. Flat string regions will therefore 
be formed, bounded by the whole of quark or anti-quark momenta, but half of gluon 
momenta. Since hadrons form from this field, the phase space is calculated by adding up 
the lengths of typical hyperbolae in each of these flat string regions. For instance, the 
invariant mass of the region involving the quark and the gluon, for the system in Eq. (|7|), 
is (fc qU ark + l^giuon) 2 = , and its contribution to the phase space for hadronisation is 



SQl 

2m, 





Such a factor was not necessary in Eq. (^). A gluon is shared between two dipoles, 
but both of those dipoles are modified by the emission of a single gluon from either of 
these dipoles. This is because the gluon in common between them would receive a recoil 
and that would change the invariant masses of both the dipoles. Formation of a single 
hadron in string fragmentation does not give recoils to the partons making up the string. 
Hadron formation occurs from the string, which has a spatial extent. Gluon emission, on 
the other hand, is a perturbative phenomenon and the new emission is thought to come 
directly from one or the other of the two partons making a dipole, or from the dipole as a 
whole regarded as a unit. Therefore, in the implementation of the dipole cascade model 
in Ariadne, the full gluon momentum is used for both the dipoles involving the gluon. 
There is no problem in conserving energy-momentum because both the dipoles can not 
emit "simultaneously" . As soon as there is an emission from one dipole, the whole system 
is updated into a new chain of dipoles. 

The expression in Eq. ([I]) does not make sense for soft or collinear gluons, as these 
would give negative contributions to the effective rapidity range. An infrared stable 
generalisation applicable for arbitrarily many gluons was presented in || and its relation 



to the hadronisation process was discussed in | I3| . For a state with a single gluon we 
replace the expression in Eq. (|7|) with: 



y m\ AtUq 
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This is a nicely interpolating expression between the results for a string with and without a 
gluon. It can be generalised to an arbitrary number of gluons by introducing the following 
quantities, defined as integrals along the directrix A.(£) (with £ a suitable parameter, here 
taken to be the energy along the directrix): 



Etot r£i 

d£i / d& 
Jo 



dA dA 

6j-i 



€2 r£,z 

d£ 3 / d£ 4 
jo 

dA dA 



dA dA 

df 3 d£ 4 



(9) 



This corresponds to a 2j-fold partitioning of the directrix into non-overlapping pieces and 
then to multiplying the adjacent vector differentials two by two. For a massless relativistic 
string, a point on the string (parametrised by the energy (a) between the point and the 
quark end) at time t will be at x(a,t): 



x(a, t) 



[A(t + (x) + A{t - a)] 



(10) 



i.e. it is given by one left-moving and one right-moving vector defined by points on the 
directrix. As the differentials of the directrix are lightlike the quantities dA(£) • dA(£') 
occurring in the integrals in Eq. (|9|) are surface elements on the string. The integrals 
in this way correspond to all the possibilities of obtaining j-fold partitionings into non- 
overlapping areas on the string surface. 

We may then define a functional T = exp(A) as 



t 



00 

1+ S W 



i.e. as the generating function of the partitioning functionals tj in Eq. (|9|) with m an 
energy scale to make it dimensionless. We note that the integrals tj will vanish if j > n+1 
where n is the number of gluons emitted. The term involving one power of m is always 
equal to and the last (nonvanishing) term will have the generic form 



1 %L S12 S23 

2 2ml ^ m o ^ m o 2mo 



(12) 



with + 1) indices of colour-consecutive partons ( index for the q, indices 1 .. n for 
gluons, and n+ 1 for the q). Besides the addition of 1 and a factor of 1/2 (both introduced 
to obtain simple factorisation properties, cf. below) the T functional in Eq. (|ll|) coincides 
with the argument of the logarithm in Eq. (^). 

It is useful to define a generalisation, T(£), by exchanging the upper limit in Eq. (|9]) 
from E tot to £ P| . This functional fulfils 



JO 



and together with the four-vector valued functional qj^. 

1 rt 



T(0 ^0 



dA(£)T(0 



(13) 



(14) 



5 Called q in § 
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we obtain the differential equations 



dT = 



TTLq 

dq T = dA 5— qx 

dql = 2(1 -%)q T dA (15) 



They can be integrated to 



T(f) = exp 



10 UIq 

,2 , 



Finally, we define a four-vector valued function X a = A u — qrjd so that 

(q T dA) 

dX = k — q T (17) 

i.e. the vector qT is the tangent to the A^curve reaching out to the directrix. From Eq. 



(|16|) we conclude that the functional T is the exponential of an area (note that the vector 
dA is lightlike, and qxdA is an area element between the directrix and the A-curve). This 
area is proportional to the functional A as mentioned in the Introduction. We also note 
that the vector qx quickly reaches the length mo- Integrating the differential equations for 
a directrix built up from a set of lightlike vectors {kj} we obtain the iterative equations 

@: 

T 

Tj+i = ( 1 + QT j k j+ i/ml)T j = 1 



7i+i 



q Tj+1 = l j+ iq T j + 1 + 2 /j+l fc i+i ( 18 ) 

The starting values are T = 1 and g^o — k . 

The factorisability of the functional T means that its logarithm, A, can be written as 
a sum containing one term for each kf 

n+1 n+1 

A = E AA . = E ln i 1 + qTj-ikj/ml) (19) 

3=1 3=1 

The AAj defined in this way is the size of a subarea just as the total A represents the 
area spanned between the X and the directrix curves. In Fig. [I] we exhibit such a region 
bordered by the "initial" qrj-i and the "final" qxj together with a hyperbolic segment 
from the A'-curve and the gluon energy-momentum vector kj from the directrix. 

If we define 70 = 1 and the vectors kj = \kj{l + 7j), then it is easy to see that the 
vector qTj is a weighted sum of the ki vectors (/ < j). The weight for any vector k m is 
exponentially suppressed with J2i= m A-^+i- 

To see that AA^ as defined in Eq. (|19D actually has the meaning of a rapidity we go 
to the rest frame of the vector qrj-i and assume for simplicity that it has the length mo- 

7 Called X in § 



7 



Figure 1: A plaquette bordered by qrj, qrj+i , kj and a hyperbolic segment from the 
<Y-curve, as described in the text. 



Further assume that the lightlike vector kj is directed along the 3-axis with the length u 
so that AXj = ln(l + uj/ttiq). Then it is easy to show that the "next" qr- vector q?j will 
have the lightcone coordinates along the 3-axis equal to (m exp(AAj), m exp(— AA 3 -)), 
i.e. the qr vector has been accelerated from rest to the rapidity AXj by "the step" kj 
along the directrix. The corresponding segment of the A'-curve is a part of a hyperbola. 

For sufficiently well ordered emissions and a resolution scale much smaller than the 
smallest transverse momentum present, the last term (Eq. (0)) will dominate in Eq. 
(PI) and therefore A(m = \k±_) will be approximately equal to £(k±) defined in Eq. (|5|). 

From the definition in Eq. fll~9|) it is easy to show (cf. the definition of the vectors qxj 
in Eq. (|T8"D) that the argument in the logarithm of the AAj is equal to (1 + Sj-ij/Arnfy 
plus terms down by at least a factor of 7. From the definition of the local k t variables, 
Eq. (0), we may also conclude that 



k% < s J^l£M+l (20) 
K s j-i,j + S j,j+V 



and consequently we obtain 



ArriQ 



k 2 



ln(-f) < n,i.,(/, ;,./,.,.;) (21) 

K ± 

This implies that there is a minimum positive value of AA, if a procedure to exclude or 
filter out emissions creating gluons with kt < k± has been used. 



3 The Properties of the AA Distributions 

We will now exhibit a very noticeable regularity of the dipole cascades, cf. Fig. |2]. In this 
figure, we show a plot of the average number of gluons as a function of k±, n(k±), versus 
£(k±) for e + e~-events with yfs = 200 GeV and for a range of values of k± from 4 down 
to 1 GeV. We have switched off the gluon splitting process so that we end up with one 
colour connected set of partons. 

There are several interesting properties of this plot: 
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• There is a linear correlation between the average gluon multiplicity (n(k±)) and the 
phase space variable £(k±) inside a very narrow band. The fact that the different 
lines obtained for different values of k± are so close to each other means that the 
slope of this linear correlation is rather insensitive to changes in the ordering variable 
in this range. It should be noted that the £(fc^)-distributions are widely varying for 
these fcj_-values as we show in Fig. |3|. We have used a running coupling (with 
Aqcd = 0.22 GeV, i.e. the default value in Ariadne). The region with £(k±) > 20 
is mainly populated by the results from small fc^-values, k± < 2 GeV. 

• We find that the linear correlation in Fig. |2| is independent of the cms energy (we 
have used different cms energies ranging from 40 GeV and up to even 2000 GeV, 
though for small energies there are too few gluons emitted to allow detailed studies), 
cf. Fig. |j. Not only do we find a linear correlation at each value of cms energy 
(when the number of gluons is more than 2), the slope of the correlation varies 
rather weakly with respect to energy. This property is also not a function of the 
total thrust, or the presence of particularly hard emissions (e.g. gluon of k± > 7 
GeV for = 200 GeV). There is a minimum value of £ ~ 10, corresponding to 
no gluon emission (for k± — 1 GeV and y/s = 200 GeV), and there is a (small) 
transition region for a few more units, until the number of gluons exceeds two. But 
after that dn/d£ is a constant. 

• As a consistency check we have used two procedures for sampling the events. It 
is possible to run the parton cascade in Ariadne down to a certain value of the 
ordering variable k±, i.e. emit all gluons with larger fcj_'s, sample the emerging 
directrix and investigate its properties, e.g. calculate the number of gluons emitted 
to that point, n(k±), the total length of the X-cmve, £(k±), the individual AA(m )'s 
and the local k t J s. After that we may continue the cascade for the same event down 
to smaller /c^-values and do the same exercise. Another possibility is to run the 
cascade down to a certain value k± = k± cut and collect a large sample of events to 
investigate. After this we can pick a new value of k± cut and repeat the process. The 
result in Fig. |2| does not depend upon which of the two methods for sampling the 
events we use. 

We observe the same effect in the parton shower in Pythia as in Ariadne. The 
ordering variable in Pythia is the parton virtuality Q 2 , and the variable k\ = z(l — z)Q 2 
(where z is the lightcone fraction taken by the emitted parton) is used as the argument 
for the running coupling. This is a good approximation for the transverse momentum in 
the case of massless parton kinematics. However, since transverse momentum itself is not 
the ordering variable in this case, there is no real lower limit to it. Therefore, in order to 
make the comparisons we filtered the events obtained from Pythia and combined partons 
such that all gluons had k t above a certain value proportional to Q. 

It is then possible to use Pythia to obtain an inclusive sample of n g (CQ) and £(CQ). 
We just note that it is possible to tune the constant C such that, as shown in Fig. |^, 
there is the same dn/d£ as in Ariadne. We find that C ~ 0.5 

It is possible to switch off the coherence conditions in the parton shower in Pythia, 
since the strong angular ordering is implemented only after the generation using the 
ordering variable Q 2 . We show the effect of this feature in Fig. [5[ For the same C we 
obtain many more partons per unit phase space for an incoherent shower, although they 
still show a linear correlation to the ^-measure. 
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Figure 2: The figure shows a plot of (n(k±)) versus l(kj_) for an e + e -event with 
200 GeV from Ariadne. 
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Figure 3: The £(/cj_)-distribution for different /cj_-values at y/s = 200 GeV. 
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Figure 4: The figure shows a plot of (n(k±)) versus £(k±) (k± = 1.0 GeV) for different 
cms energies from Ariadne. 
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Figure 5: The figure to the left shows a plot of (n(k±)) versus l(kj_) for an e + e~- 
event with ^/s = 200 GeV from Pythia. The figure to the right shows (n(k±)) ± a 
versus £(k±) at ^Js = 200 GeV, k± = 1.0 GeV for coherent and incoherent showers 
in Pythia. 
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Figure 6: Scatter plot of two adjacent A£ values 



The linear relation between (n) and £ can be understood as follows. If one considers 
a variable such as £ defined as in Eq. (|5|) where the parts A£j are distributed with the 
averages fij and the variances crj, then £ itself will mostly have a gaussian distribution 
according to the central limit theorem. The average is p — and the variance a 2 = 
Y^ffj- The numbers p and a 2 are in general well-defined for large n unless the values of 
the A£j are correlated. We find that the Aij's indeed show very little correlation, even 
between adjacent dipoles, cf. Fig ||. We also find that the average \ij shows only a weak 
dependence upon the total number of gluons in the event, n, for a wide range of n. 

We will therefore go over to the investigation of the distribution of the A£j defined 
in Eq. (H). We notice that the first and the last term in Eq. (|j) are different from the 
other terms in one important aspect. They are contributions from dipoles involving a 
quark or an antiquark (with different colour charges than gluons), whereas all the other 
terms come from purely gluonic dipoles. It turns out that the dipoles at the ends have 
slightly different distributions than all the rest. However, if the gluon splitting process is 
switched off, there are only two such dipoles in an event. This means that we can write 
for the average p = 2(/z q ) + {n — l)(fi gg ). In Fig. | we see this effect clearly, because it 
is only for n > 2 that we find a constant slope, i.e. when we have purely gluonic dipoles 
present. The properties we noted above are about variation of the number of gluons or 
dipoles with respect to £(k±_) and the stability of the slope of the correlations, and should 
correspond to the properties of purely gluonic dipoles. Therefore, from now on we will not 
include the very first and the very last A£-values in the investigations. We have seen that 
purely gluonic dipoles, irrespective of their position along the directrix, have the same 
distribution in size. In Fig. [7|, we show the distribution = Pi(A£) from Ariadne 
and observe that: 

• The distribution Pi has an average \x and a variance a 2 , consistent with the results 
in Fig. |2|. We have found that it is independent of the cms energy, of the global 
event variables like thrust and the total £- value for the event, and also of hard gluon 
emissions. 
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Figure 7: The = ^^-distribution from ARIADNE. 



• The distribution Pi is rather insensitive to the value of k±. In Fig. ^ we show the 
average value of A£ as a function of k = ln(fc^_/A 2 ) and we find that it changes 
about 5 %, when the ordering variable goes from 0.8 to 3 GeV. 

• There are no noticeable correlations between adjacent Ats as we have noted before. 
In Fig. ^, we show the scatter plot beween two adjacent A£, taken from stochas- 
tically chosen pairs in many different events. In Fig. |9|, we show the values of the 
average a{n) (defined as a(n) = (J2]=iA£j)), and the standard deviations around 
it, for chains of n connected A£. For this figure the statistics was collected from 
events with a fixed number of gluons. 

• In order to further investigate the independence of the A£ values in an event we 
have examined events containing a fixed number (JV) of gluons. We may then ask 
about the multiplicity distribution, P(Na), of events with values of A£ < A. 
We find a binomial distribution with the mean Np and variance Np(l — p) where 
p = / A dxPe(x) as expected from uncorrelated dipoles. We have also investigated 
the average number of A£- values, which are consecutive and satisfy A£ < A. We 
find a geometric distribution, weighted with the average number of A£ not fulfilling 
the condition (which is, once again, expected when we have uncorrelated dipoles). 
We also find that the distribution P^ n obtained from the sum of the lengths of n 
consecutive A£ is an ra-fold convolution of the distribution P^. 

The corresponding results in Pythia are very similar to the results in Ariadne (with 
the appropriate choice of the constant C as discussed above). To the left in Fig. [l(| , we 
show the distributions from Pythia and note that Pythia is somewhat more concen- 
trated towards small A£ values and also has a somewhat wider tail so that the average 
values are the same as in Ariadne. The "incoherent" option gives a much more narrow 
distribution centred around smaller values, as expected. 
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K 

Figure 8: The average value of A£ as a function of k = ln(/c5_/A 2 ). Note that the variation 
in the average A£ in the figure is much smaller than the width of the distribution of A£ 
as seen in Fig. 0. 



It is also possible to investigate the distributions in terms of the A-measure defined 
in Eq. (P"§1). In Fig. |TD| to the right, we show two distributions (one from ARIADNE 
and one from Pythia), to be called P\(AX) using as before 2mo = k± = CQ where mo 
is the argument of the A-measure. Once again we find the same properties as for the 
distributions Pe(A£). 

We conclude that, inside the whole region relevant to the hadronisation procedure in 
String Fragmentation, the partonic states are dominated by this structure of independent 
entities. It is tempting to consider them as a kind of collective variables for the QCD force 
fields. The A£ variables represent the local size of the dipoles between adjacent parton 
energy-momentum vectors, whereas the AA values correspond to a similar although more 
non-local property depending upon a set of adjacent vectors. 

These properties of the distributions of Ats are properties of gluonic dipoles, and 
result from perturbative QCD. In the next section we will make models in the LLA and 
the MLLA schemes based on simple gain loss considerations on the dipole cascade model, 
and show that it is possible to understand many qualitative features of these distributions. 

The AA's on the other hand are complicated objects, and there is no obvious way to 
calculate their distributions except by relying on the similarity between the A-measure 
and £. But the fact that there are as many AA's as Ai'a, and the AA's have parallels 
for each of the properties of the A£'s discussed here, lead us to introduce the notion of 
Generalised Dipoles (GD) to be the "source" of the contributions to the AA/s. Whereas 
a dipole is spanned between two colour connected partons, a GD corresponds to one flat 
block or "plaquette" in the surface traced between the directrix and the A-curve by the 
tangent to the A'-curve. Regularities in the properties of dipoles which are mirrored in 
the GD's would be passed on to the hadronised state in a hadronisation scheme based on 
the A-measure, such as the one introduced in |J. 
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Figure 9: The figure shows the average A£ (and the standard deviations around it) for 
chains of n connected A£ called a(n) (a(n) = (S"=i ^£j)). The statistics is collected from 
events with a fixed number of gluons. 
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Figure 10: The figure to the left shows the P^(A£)-distribution from Pythia com- 
pared to the distribution obtained from Ariadne. The figure to the right shows 
the corresponding plot for the P\(A\) distribution. 



4 Models to Describe The Dipole Distributions 



4.1 A Simple Cascade Model 



Multiplicity distributions were presented in the parton cascade formalism in [14| and in 
the dipole cascade formalism in |Tj| to leading logarithmic order. In this paper we will 
follow the approach presented in [TE\. 

In order to provide a theoretical framework for our findingsin the previous section, we 
will consider an analytical model for a dipole cascade. We will assume that there is a 
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distribution, /(//, K)d[i, that describes the dipoles with size \i at the scale k = ln(fc^/A 2 ) 
(where k± is the ordering variable and u = \n(M^/kj_), in terms of the dipole masses 
Mrf). We will consider the change in / when we make a small change in the scale k. We 

note that jj + k = ln^4- is independent of k. If dipoles never decayed, the distribution 
/(//, k) would not change if \i + k is kept constant as k varies. Therefore any change in 
the distribution along the lines fi + k = constant must purely come from decays. 

There are two contributions to this. A dipole of size /i may decay into two smaller 
dipoles. A dipole of size bigger than /i may decay such that one of the daughter dipoles 
is of size \i. Gathering these terms together we get a partial integro-differential equation, 
cf. [[Tj], [T7[] , based upon Eq. [I] with a. e s = for gluonic dipoles (cto = u _ 2 N f /N ) : 



l*f - 2 / dfJf(n', 

J a 



[22) 



The factor of 2 in the contribution from the larger dipoles stems from the fact that we are 
only considering "central" or purely gluonic dipoles here. When a gluonic dipole decays, 
there are 2 equivalent ways to obtain a gluonic dipole of size \i. 

This integro-differential equation can be made into two coupled differential equations 
in terms of the first two moments of /, Ni and 

Nj(fjL,K) = / tfy4i(f4i) J_1 /0*Lj K ) 

DN l = - — (N 2 -2aN l ) 

DN 2 = -A^ + — a 2 N x (23) 

K 

For the normalisation we note that Ni(0,k) is the total number of dipoles available at 
the scale k, to be called n in this section. A^ 2 (0, k) is their combined length, to be called 
I. The equations for these quantities can be obtained directly from Eqs. ( p3|) in the 
LLA-scheme: 

dn a - 
dn K 

Tk = ~ n (24) 

These equations are solvable in terms of combinations of the modified Bessel functions 
(cf. the next subsection) as soon as we have defined a "starting value", n max = L , for 
the cascade. 

In order to solve the coupled equations in Eq. (|23|), we introduce the combination 
g(p,, k) = N 2 — fiNi, which fulfils the equation: 

Dg = ^(w) (25) 

K 

with the boundary value 



It is related to Ni through Ni(p, k) = —i-g(fji, k). 



g(0,x) = £(x) (26) 
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Including the boundary condition above, we obtain: 



d_ 

dfj, 



£{jx + k) exp 



dy—(fi + K - y) 

y 



(27) 



This is a derivative of two contributing factors. The first is I evaluated at // + k = 
ln(Mj/A 2 ), i.e. at the largest k value where the dipole could have been produced. It is 
multiplied by the probability that the dipole has not decayed until k. In field theoretical 
language this second contribution, the Sudakov form factor, corresponds to the sum of all 
the virtual corrections during the "lifetime" of the dipole. This Sudakov form factor can 
be easily calculated 



exp 



/ dy—({i + K-y) 
Jk y 



exp(a /x) 
'1 _|_ H)<*o (/•»+«) 



It fulfils 



5(0, k) 
dS 



:28i 



(29) 



-a \n(l + ^)S 

K 



The derivative vanishes at fi = 0. 

The distribution / is obtained by a partial differentiation of Ni, cf. Eq. (EJ 



so that using Eqs. 



dfx 

and (^) we obtain: 

= [£(L Q ,n + K)a \n(l + fi/n) 2 + 

2n(L , n + K)a ln(l + /// k)]S((/,, k) 



(30) 



(31) 



The result can be understood in a simple and useful way. The logarithmic factors in 
the parenthesis are equal to 



a ln(l + /jl/k) = 

J K 



m+ k a dx 



x 



(32) 



i.e. the probability that there will be a dipole breakup somewhere between /i + k, the 
largest virtuality where a dipole with size \i at k can be produced, and k, where it is 
found. There are two such factors multiplying all interior points in the dipoles at \i + k 
, counted by £((/, + k), and there is one factor multiplying twice the number of dipoles at 
fi + k (because if there is already a dipole one can keep either its "left" side or its "right" 
side and produce the "other" side of the dipole with size \i by a breakup in between). 
Finally, the Sudakov form factor describes the probability that there is no decay affecting 
the existence of the dipole at k. 

In the limit when Lq ^> (// + k), and /i + k is significantly positive, the two functions 
n and I factorise and we may write in the LLA 



n(L ,x) ~ V(Lq)(x) 4 exp(— 2y/a x) 
£(L ,x) ~ Jx/a n(L ,x) 



(33) 
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Figure 11: The function, obtained from Eq. |31], normalised to unity. 



In Fig. |TT], we show the function f K (fi), obtained from / in Eq. (j3lD by normalising it 
to unity, for a set of values of the parameter k. There is a weak dependence upon the 
ordering variable k and for large /^-values it has an essentially exponential falloff with a 
slope very similar to the ones we find in the distributions Pe(A£) and P\(A\) in Section 
i 

For small values of \i there are quantitative differences. The reason is that the LLA 
is not sufficient to describe the behaviour of small dipoles. In that case it is necessary to 
include both the contributions from the {J2 P ) term in the dipole cross section in Eq. ([]]), 
and also to take into account the recoil contributions which are particularly noticeable 
for small dipoles. We will consider this question in the next subsection. 

4.2 The MLLA and Earlier Investigations 

In the so-called Modified Leading Log Approximation, subleading effects from the polaris- 
taion sum are included ]7|, [TJ| . In general, the assumptions have been that the cms energy 



y/s (or rather its logarithm L = ln(s/A 2 )) is very large and that the main contributions 
from the cascades will come from well-ordered emissions, i.e. that k±j ^> k±j + \ (at least 
for adjacent emissions). 

Using this approximation, it is possible to write out analytical equations for the ex- 
pected changes in the average i (to be called £(k)) and the average n (similarly n{n)) for 
a given ordering variable k = \n.{k\/ 'A 2 ) [0, [17|, |18| : 

M_ 

fjf) _ 

— = - aeS (£-nS) (34) 

UK 

Eqs. ( |34| ) are simplified versions of the Modified Leading Log Approximation (MLLA) 
formulae. A more elaborate treatment should include the difference between the quark and 
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antiquark endpoint dipoles with a different effective coupling a and a different correction 
5, cf. ||18|| , but these are small corrections if the starting value L = ln(s/A 2 ) is large. 

The interpretation is that in a step dn, 1 changes by n contributions from the scale 
change according to Eq. (|5|), just like the second line in Eq. (p4|) . Further, the change in 
n in that step is given by the coupling times the available phase space. This phase space is 
given by i to first approximation; but there is a correction from each dipole from the term 
J2 P m Eq. ([!]). The polarisation sum J2 P — 2 near the centre of a dipole but decreases 



in the neighbourhood of the emitters. It has been shown |18| that the suppression of 
emissions because of this can be approximated by (and this is essentially the MLLA) a 
constant decrease in the dipole size 5 = 11/6 for gluonic dipoles (this number only depends 
upon the triple gauge boson vertex for 3 colours). In Eqs. (|34D it is also assumed that 
the gluon splitting process g — > qq is neglected and therefore a c s = oco/k with «o = 6/11 
for all dipoles and not only the purely gluonic ones. 

The decay of one dipole also affects the two adjacent ones due to recoils. A method 



to estimate this effect was suggested in ||19|| . Therefore an extra "loss-term" is added to 
the first line in Eq. (|34|): 

dl _ dn 

— = -n- C r — (35) 
an an 

with a constant recoil correction estimated to be C r ~ 2. 

These equations can be solved using the following combinations of the modified Bessel 
functions: 

l x {x) = v / 2x 7/2 J 7 (2 V / o^) 

I 2 (x) = v / 2^a; (7 - 1)/2 V 1 (2 v ^) (36) 

and similarly for /Ci and K, 2 in terms of the exponentially falling Bessel function K. Note 
that they are normalised so that 

X!(x)/C 2 (x) +l 2 (x))C 1 (x) = x 1 - 1 (37) 

and that there are simple differential relations between them: 



dx 

dI 2 a (7-I) 

— — = — Xi(x)H M{x) (38) 

dx x x 

Similar relations hold for the exponentially falling modified Bessel function pair (/Ci, — A^)- 
If we start the cascade at a large value of L = ln(s/A 2 ) > k we obtain, with the 
boundary values i = and n — 1 at L, the results: 

a = L 1 ^ [C r IC 2 (L) -JCx(L)} (39) 

b = L^lCM^+l^L)} 

n = 6/C 2 (k) — aX 2 {K) 

I + C r n = ali (k) + bid (k) 

using the notations 7 = 1 + ao(5 + C T ). 

It is difficult to obtain an anlytically solvable partial integro-differential equation corre- 
sponding to Eq. (^2|). The method of moments runs into difficulties because the variation 
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in the J2p term in the cross section (Eq. ([!])) affects the integral contribution, i.e. the 
gain from the decays of larger dipoles, in such a way that the equations for the first few 
moments are no longer independent of higher moments. Therefore, we will directly go 
over to the solution in Eq. (|3l|) instead, and modify the ingoing terms according to the 
MLLA, keeping in mind the interpretation given there. 

We assume that each dipole of size /i can decay only in the interior excluding a region 
5/2 on each extreme. Then the effective phase space at (/x + k) in Eq. (|3l|) is changed 
from £(fjb + k) into [£(p + k) — 5 ■ n(/x + n)] . A dipole of size \i is produced by a breakup first 
at a point k — K\ and then at k 2 < K\. The region where the "no emission probability" 
is to be calculated in the Sudakov form factor is changed in the MLLA, first by a loss of 
5/2 at K\ and then by a loss of a further 5/2 at k 2 . 

( f K2 aodz(fi + k — 5 — z) \ I r^i a dz(fi + k — 5/2 — z)\ 

5ki ' K2 = exp \i z J exp v A 2 I ) x 

6 



exp 



J k\ Z I \ K / 



where S(fi, k) is the Sudakov form factor without the term 5 as given in Eq. fl28|). Since 
this is a function of K\ and k 2 now, we obtain, for the contribution to /(/J, k) from 
situations where both gluons involved in a dipole are emitted below /j + k, the result: 

/i = 2(£ - 5n)(j2 + k) / S KuK2 



4 ((1+ -I)' 

= (E-6n)(ji + K)-± *-S(ji,k) (41) 

5 Z 

The fact that the first break may be either on the "left" or on the "right" side contributes 
a factor of 2 in the above equation . 

In the same way, for the second contribution related to the number of dipoles at fi + K, 
ie the situations where only one gluon involved in a dipole is produced below the scale 
ji + k, we obtain, 



4 (1 + (m-*/2)/k)-*--1 

h = 0(ti - 5/2)n(jJi + k)-^ '-Sin - 5/2, «) (42) 

5 

In this case there can never be a breakup until the dipole size has at least reached 5/2, 
because one of the end points is already fixed at the virtuality /j + k. 

We note for consistency that the result f s = f 1 + f 2 has the limit / in Eq. (|3lD 



when 5 vanishes. In Fig. [12], we show the result of this MLLA version of the dipole size 
distribution. We have used a fixed value 5 = 11/6 but we have allowed the parameters a 
and C r vary. The reason for allowing «o to vary is that it is feasible to shut off the gluon 
splitting process in Ariadne, but the running coupling will always receive contributions 
from the different number of "active" flavours at different virtualities. From the figure to 
the left, we conclude that it is in between the distributions from Ariadne and Pythia 
for the central dipoles when we use a = 0.70 (i.e. for three to four active flavours), and 
C r = 2 (according to |T9|]). This result is rather insensitive to C r and it can therefore not 



be used as a test for the claims in |]T9|. We have also checked the dependence upon the 
ordering variable. In Fig. to the right, we show that f$ also exhibits a slow change of 
its average value as a function of the (relevant) ordering variable like the Ariadne and 
Pythia P f distributions. 
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Figure 12: The figure to the left shows the distribution /#(//,/«), obtained from 
the MLLA approach at k± —IGeV, compared with the Pi(A£) distribution from 
Ariadne and Pythia (cf. Fig. |H]). The figure to the right shows the average 
value of /i as a function of k, from the same distribution fs{fi, k)- In the figure to 
the right, we also compare with the average value of i (cf. Eq. ([H])) as a function of 
k obtained from Ariadne and Pythia, and we also include the standard deviation 
in the result from Pythia and Ariadne. This result is obtained using an = 0.7 
and C r = 2 in the MLLA distribution as explained in the text. 



5 Concluding Remarks 

The number of gluons emitted above a certain value of the ordering variable shows an 
interesting linear correlation (with gaussian fluctuations) with the phase space available 
at that scale. We have called the phase space variable the ^-measure in this paper to 
distinguish it from the A-measure, which is an infrared stable generalisation of the £- 
measure, useful in the context of hadronisation. 

Both these quantities are functions of not only the resolution scale (kj_ cut for a dipole 
cascade, r in string fragmentation), but also of the precise geometry of the state under 
consideration. Therefore they change with every new gluon emitted in the system. 

The resolution parameter relevant for hadronisation is a parameter fixed by the break- 
up properties of the string and the mass spectrum of the hadrons produced. The A- 
measure is constructed in such a way that addition of gluons at transverse momenta 
much smaller than the scale T does not change its value. In that way it is infrared 
stable. 

However, since the resolution parameter for the ^-measure is the transverse momentum 
of the gluon to be emitted, the transverse momenta of all gluons already emitted must 
be higher in an ordered dipole cascade. Therefore, the only kind of infrared stability 
relevant for this measure would be the stability, when the ordering variable, which is also 
the resolution parameter, is small. Technically this is not the case. However, since in this 
approach we will never compute phase space available at a value of the ordering variable 
below the Aqcd, there is an effective infrared stability for the phase space, if all local 
transverse momenta k t j are required to be above k±. It is easily shown that the mass of 
any dipole with a gluon at one end is greater than the transverse momentum k t of that 
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gluon, which in turn, will be greater than Aqcd- Besides, we must remember that it is 
kinematically not possible for dipoles of mass smaller than 2 Aqcd to emit gluons above 
Aqcd- 

The two measures discussed here are quite strongly related. If one keeps only the last 
(often the dominant) term in the argument of the logarithm in the A-measure, one gets 
an expression very similar to the I- measure. 

The linear correlation between the ^-measure and the number of gluons is a direct 
consequence of the following facts that we have demonstrated. The ^-measure is the sum 
of a series of terms A£j, one for each dipole. A given A£j is the contribution from one 
particular dipole to the ^-measure, and is in some sense the "size" of that dipole. There 
are two kinds of dipoles. The first kind involves either the quark or the antiquark energy- 
momentum. The second kind involves only the gluon energy-momenta. These two kinds 
of dipoles have different distributions in sizes. 

It turns out that the size of the dipoles do not show any significant correlations. 
Therefore the ^-measure of a state consisting of a certain number of dipoles, which is just 
the sum of the corresponding A^-'s, is distributed like the sum of several independent 
random variables, i.e. like a gaussian with a mean that is the sum of the means of the 
individual independent distributions. This shows up as a linear correlation between the 
number of partons and the ^-measure. 

More interesting than the linear relation is the stability of the slope of these lines, 
which represents the number of gluons per unit phase space, with respect to changes in 
global event parameters like cms energy or thrust and even with changes in the ordering 
variable. It seems that parton cascades tend to form structures at one (relative) size when 
sizes are measured with the ^-measure with a resolution parameter proportional to the 
ordering variable. As we go down in the ordering variable, the running coupling should 
increase the number of emissions per unit phase space. However, the phase space variable 
itself scales appropriately to effectively absorb any increase of emissions per unit phase 
space. 

We have also made another observation in this paper. Just like the ^-measure, the 
A-measure can be thought of as a sum of terms AAj, one for each vector along the directrix 
in colour order. Out of the two vectors involved in a particular AAj, only one is a vector 
along the directrix. The other vector, in a sense a cumulative variable, is a weighted 
average of the vectors k^s along the directrix with % < j. 

It is interesting that all the properties of the ^-measure we have discussed here are also 
reflected in the A-measure. There is a similar linear correlation to the number of gluons, 
similar lack of correlation among adjacent terms and similar stability with respect to 
change of global event variables and the ordering variable. Partitioning the A-measure into 
a series of AAj's corresponds to partitioning the surface spanned between the directrix and 
the X-curve into a series of "plaquettes" or flat regions along the directrix. Each plaquette 
is spanned between a parton energy momentum vector kj, a vector qrj-i (weighted average 
of fcj's along the directrix with i < j), another vector qxj and a hyperbolic section of the 
X-curve. The size of the plaquette AAj is the length of the section of the X-curve, and is 
determined by kj and qrj-i- Since these plaquettes share so many properties with dipoles, 
we will call them "Generalised Dipoles" . 

We note in passing that the possibility of partitioning the contributions to the A- 
measure as contributions from connected flat regions has an interesting consequence for 
the fragmentation of a string according to the Lund model area law. Such a fragmentation 
scheme was presented in M, which was closely related to the A-measure and the X- 
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curve. Particle production in this scheme can be thought of as partitioning the plaquettes 
mentioned here into smaller plaquettes, one for each hadron. Since the hadron energy- 
momenta will mostly come from these flat regions in the string, particles stemming from 
one GD will have energy-momenta aligned in a plane in Minkowski space, up to transverse 
momentum fluctuations. We have found such chains of particles and examined their 
properties. These chains, which we will call "coherence chains" are important for the 
study of Bose Einstein correlations in the Lund Model for multigluon systems, which we 
will elaborate on in a forthcoming publication. 

We have also presented here an analytical model for the observations we have made 
about dipoles, in the LLA and MLLA schemes. There are limitations to this approach, 
since complications such as the polarisation sum make it very difficult to obtain differential 
equations which can be solved analytically. Nevertheless, our analysis based on simple 
gain-loss considerations broadly reproduces the qualitative features of the observations. 
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